Correlation effects in insulating surface nanostructures 
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We study the role of static and dynamical Coulomb correlation effects on the electronic and 
magnetic properties of individual Mn, Fe and Co adatoms deposited on the CuN surface. For 
these purposes, we construct a realistic Anderson model, solve it by using finite-temperature exact 
diagonalization method and compare the calculated one-particle spectral functions with the LDA+t/ 
densities of states. In contrast to Mn/CuN and Fe/CuN, the cobalt system tends to form the 
electronic excitations at the Fermi level. Based on the calculated magnetic response functions, the 
relative relaxation times for the magnetic moments of impurity orbit als are estimated. To study the 
effect of the dynamical correlations on the exchange interaction in nanoclusters, we solve the two- 
impurity Anderson model for the Mn dimer on the CuN surface. It is found that the experimental 
exchange interaction can be well reproduced by employing [7=3 eV, which is two times smaller 
than the value used in static mean- field LDA+[/ calculations. This suggests on important role of 
dynamical correlations in the interaction between adatoms on a surface. 

PACS numbers: 71.27.+a, 73.20.At, 75.20.Hr 



Recent scanning tunneling microscopy (STM) experi- 
ments on atoms deposited on surfaces have demonstrated 
many interesting and promising results varying from fun- 
damental phenomena, such as the Kondo effect^ or gi- 
ant magnetic anisotropy,^ to practical realization of logic 
gates. ^ Theoretical support of these experiments requires 
the construction and solution of realistic models that can 
capture the most essential couplings and effects in the 
system. In case of the impurity deposited on an insu- 
lating surface, such as Mn (Co, Fe) adatom on the CuN 
surface, it is natural to start the study with a spin-type 
Hamiltonian^i^ where the impurity is considered as a lo- 
calized spin. This simple and attractive approach has 
been successfully applied to reproduce the Kondo effect 
in Co/CuN^ or step-type excitations in the conductance 
spectra observed for Mn/CuN."^'^ 

The next step toward realistic modeling of Mn/CuN, 
Fe/CuN and Co/CuN systems would be to take into ac- 
count the orbital nature of the adatom and hybridiza- 
tion effects between the adatom and surface states. This 
can be done by using a density functional theory (DFT) 
method^ or a many-body Anderson impurity model ap- 
proach, or their combination. The first-principles DFT- 
based calculations provide an important information con- 
cerning ground state properties of the surface nanosys- 
tems. For instance, the densities of the states obtained 
for Mn, Fe, Co/CuN systems by using the local den- 
sity approximation (LDA), generalized gradient approx- 
imation (GGA) and local density approximation taking 
into account the on-site Coulomb interaction (LDA+/7) 
demonstrated that the one-particle excitation picture de- 
pends on the symmetry of the particular Zd orbital of the 
adatom/ 

On the other hand for the simulation of the electronic 



or magnetic excitations in a surface nanosystem, the An- 
derson model approach seems to be the most natural 
choice. The parameters for this model can be estimated 
by using the first-principles calculations. However, the 
solution of the Anderson model in case of the five-orbital 
impurity at experimental temperatures is still a challeng- 
ing methodological and technical problem. 

The correct description of the exchange interactions 
between impurity magnetic moments is another compu- 
tational and methodological problem. Numerous inves- 
tigations based on the LDA+/7 or GGA+/7 approxima- 
tions clarified the role of the static Coulomb interactions. 
The information about the influence of the dynamical 
Coulomb correlations on magnetic excitations in surface 
nanostructures is still limited. For instance, the authors 
of Refjsl proposed that a disagreement between theoreti- 
cal and experimental estimates of the isotropic exchange 
interactions can be related to the lack of the dynam- 
ical correlation effects^ that renormalize the electronic 
spectrum near the Fermi level. However, the proper mi- 
croscopic mechanisms of correlated exchange process are 
still unknown. 

In this paper, we address above mentioned problems 
and study the correlation effects in transition metal 
atoms deposited on insulating CuN surface. We first an- 
alyze the spectral functions obtained from the LDA+/7 
and Anderson model calculations. Based on the calcu- 
lated spin- flip- type response functions of the impurity or- 
bitals we discuss the origin of the low-energy excitations 
and estimate the relaxation times for magnetic moments 
of impurity orbit als. We clearly observe the difference in 
magnetic excitations between Co and Fe (Mn) systems. 
In contrast to the other, the cobalt system demonstrates 
an orbital polarization of the response functions. To con- 




FIG. 1. (Color online) Schematic representation of the in- 
dividual transition metal atom (Mn, Fe and Co) on the CuN 
surface (top view) simulated in this work. 



TABLE I. Structural parameters for Mn, Fe, and Co adatoms 
supported on the CuN surface obtained from the LDA and 
LDA+[/ ([/=3, 6 eV) calculations: the surface- ad atom dis- 
tance (in A) and the angle (in degrees) of nitrogen-metal- 
nitrogen bond. See text for details. 



Mn 



Fe 



Co 



LDA 
LDA+[/ ([/=3 eV) 
LDA+[/ ([/=6 eV) 



L49 (118) 1.46 (118) 1.43 (118) 
1.46 (123) 1.38 (129) 1.37 (125) 
1.53 (118) 1.43 (126) 1.40 (123) 



nect lifetimes of quantum spin- flip excitations with relax- 
ation times of orbital magnetic moment, we used a com- 
bination of linear response theory and the Bloch equa- 
tions approach. The effect of the dynamical correlations 
on the magnetic exchange interaction is demonstrated on 
the example of Mn dimer on the CuN surface. 

Methods. In our study, we use two complementary 
approaches to investigate the electronic structure and 
magnetic properties of transition metal atoms deposited 
on the CuN surface: (i) the projected augmented-wave 
(PAW) method^^ for individual Mn, Fe, or Co atoms, 
and (ii) the tight-binding linear-muffin-tin-orbital atomic 
sphere approximation (TB-LMTO-ASA) method^^ for 
Mn dimer. The exchange and correlation effects have 
been taken into account by using the LDA and LDA+/7— 
approaches. The computational details are presented in 
the Supplementary materials. ^^ 

Geometry. The CuN(lOO) surface is built up from a 
(100) surface of fee copper and nitrogen atoms uniformly 
embedded into the topmost layer with the Cu:N ratio of 
2:1 (Fig.l). The formation of covalent bonds in the Cu-N 
network favors an insulating behavior of this surfaced 
The deposition of a single transition metal (TM) atom 
on top of a Cu atom results in a slight distortion in the 
surrounding geometry due to formation of the N-TM-N 
bonds. ^^ Equilibrium geometry for different single atom 
adsorbates can be characterized by the distance from the 
adatom to the unperturbed topmost surface layer, and by 
the N-TM-N angle. In Table IH we list these quantities 
for Mn, Fe, and Co adatoms considered in this work. As 
one can see, in all three cases the geometries are similar, 
which suggests the same bonding mechanism for different 



adatoms. In contrast to the dimer and longer chains of 
TM atoms on CuN,^^ the on-site Coulomb correlations 
do not significantly affect the geometry of the system in 
the presence of a single atom on the surface. 
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FIG. 2. (Color online) Partial one-particle spectral func- 
tions of Mn/CuN (red), Fe/CuN (gray) and Co/CuN (blue) 
obtained from LDA calculations. 

LDA results. First, we analyze spectral functions cal- 
culated by using the local density approximation (Fig. 2). 
All the systems demonstrate similar one-particle excita- 
tion spectra with the shifted chemical potential due to 
the different number of 3d electrons. The main feature 
of each spectrum is a peak at the Fermi level that orig- 
inates from the hybridization between the adatom and 
surface states. The width of the LDA peak of about 1 eV 
is approximately the same as for Co adatom on Pt(lll), 
where strong hybridization and correlation effects were 
observed. ^^ Thus the straightforward transition from the 
complex electronic Hamiltonian to a simple spin model 
is not justified. 

To give a quantitative estimate for LDA results and 
to extract the parameters for the many-body Anderson 
model, one can use the different schemes. One commonly 
used way is a projection procedure^^ that defines the 
Wannier functions and the corresponding parameters for 
a tight-binding Hamiltonian. However, the energy bands 
of the surface nanosystems are strongly entangled in a 
wide range, which makes the band structure rather com- 
plicated. Moreover, since we use the exact diagonaliza- 
tion solver there is a limitation on the number of the 
surface orbit als. It is not clear which surface orbitals 
should be included in the Wannier function (WF) basis. 



Thus the construction of a moderate-size WF basis is not 
trivial for surface nanostructures. 

A more appropriate approach is to define some effec- 
tive surface orbitals that reproduce main features of the 
LDA density of states. It can be done within a minimiza- 
tion scheme where the non-interacting Green's function 
is diagonal in the orbital space and each LDA spectrum 
is fitted by the following expression 
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where e^ (e^) is the energy of the impurity (effective sur- 
face) states and Vip is the hopping between ith impurity 
and pth effective surface orbitals. The calculated param- 
eters are presented in Ref.l3. The typical values of the 
hybridization between the surface and impurity states for 
all the systems are ranging from 80 meV to 200 meV. 

The choice of U . The correct description of the local- 
ized nature of the 3d states of a transition metal adatom 
requires the proper account of the Coulomb correlation 
effects. It can be done by using the LDA+/7 calcula- 
tion scheme. In such calculations, one needs to spec- 
ify the values of the on-site Coulomb and intra-atomic 
exchange Jh- There is a number of constrained pro- 
cedures to calculate U and Jh within the DFT-based 
methods. However, depending on the approach that can 
be local density approximation, generalized gradient ap- 
proximation or random phase approximation (RPA), the 
resulting Coulomb interaction values are rather different. 
These values are also sensitive to the details of the sur- 
face. For instance, the constrained-LDA method gives U 
= 6.75 eV for individual Co atom on Pt(lll) surface. ^^ 
On the other hand the constrained-GGA results in U=0.8 
eV in case of the Co/CuN system.^ 

We also note that a common practice is to use the U 
and Jh as fitting parameters to achieve the best agree- 
ment with known experimental data. For instance, the 
LDA-\-U method with U = 5.88 eV gives the isotropic ex- 
change interactions between Mn atoms on the CuN sur- 
face that agree well with experimental estimates. ^^ As 
we will show below, this value of the U parameter can 
be substantially renormalized by taking into account the 
dynamical Coulomb correlation effects. Thus there is no 
unique set of the U and Jh parameters for a particular 
surface nanosystem. 

In this work, we simulate the nanosurface systems in 
different interaction regimes with two sets of the U and 
Jh parameters: U = 3 eV {Jh =0.9 eV) and U = 6 
eV {Jh = 0.9 eV). The first set corresponds to nearly 
itinerant state with the effective Coulomb interaction 
{Ueff = U — Jh = 2.1 eV), whereas the second can be 
associated with a localized regime {Ueff = U — Jh = 5.1 
eV). As the U parameter of a single transition metal atom 
is known to decrease toward the surface^ two interacting 
regimes can be interpreted as two possible states, cor- 
responding to different separations of the adatom from 
the surface. Moreover, comparing the Anderson model 



results for defined U values will give us opportunity to 
retrace the corresponding changes of the electronic struc- 
ture and magnetic properties. 

Anderson model results To simulate the many-body 
correlation effects in Fe, Co and Mn adatoms on the CuN 
surface, we use the single-impurity Anderson model that 
can be written in the following form: 

per ia 

^^'^9^BBeff{ni^ - mi) 

i 
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Here e^ and e^ are the energies of the impurity and bath 
states, dl^ and c+, are the creation operators for impu- 
rity and surface electrons, respectively. B^ff is an ef- 
fective magnetic field applied in z-direction, Vi^ is the 
hopping between the impurity and surface states, Uiju 
is the Coulomb matrix element and the impurity orbital 
index i (j, /c, /) runs over the Zd states. In our work, we 
use the full form of the Coulomb interaction matrix. For 
each system under consideration the Anderson model is 
solved by using the finite-temperature exact diagonaliza- 
tion method on the real energy axis.^^ 

It is important to discuss the magnetic field contribut- 
ing to the Hamiltonian. The experimental values of the 
magnetic fields (1-7 T) being used in our model can lead 
to nearly paramagnetic state with tiny magnetization^^* 
or to magnetic solution with a low-spin state. It is mainly 
due to the use of the full Coulomb interaction matrix that 
contains spin-fiip terms destroying the magnetic state. 
To obtain the magnetic solutions with the magnetic mo- 
ment corresponding to the experimental spin, we need to 
use the effective magnetic fields of about 20-50 meV. 

To assess the effect of the dynamical correlations, we 
compare the one-particle spectral function obtained from 
the Anderson model calculations with those calculated 
with the LDA+/7 approach Fig. 3 and 4. In case of Mn 
and Fe/CuN systems, we observe the tendency of the An- 
derson model solution to enhance the splitting between 
spin- up and spin-down states for partially filled orbitals. 
It is not the case for Co/CuN system, where the account 
of the dynamical correlations for Co adatom results in 
new excitations at the Fermi level. In our model, the dif- 
ference between Mn/CuN and Fe/CuN concerns the xy 
orbital, which demonstrates a tiny magnetization in the 
case of Fe adatom. It is also important to note that we 
consider the magnetic solution of the Anderson model, 
which means that the correlation effects could be par- 
tially suppressed by magnetism. 

Spin- flip response functions. To study the magnetic 
excitations, we have calculated the spin-fiip response 
function for each impurity orbital using the Kubo's linear 



TABLE II. Comparison of the magnetic moments (in /jb) of 
Mn, Fe, Co atoms obtained from the LDA+[/ and Anderson 
model calculations with [/ = 3 eV and [/ = 6 eV. 

Exp. LDA+[/ LDA+[/ Anderson Anderson 
U=3 U=6 U=3 U=6 

Mn 5 (5'= 1)^^ 
Fe 4{S = 2)20 

Co 3{s^ ly 
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4.4 
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4.7 
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3.2 


3.7 


3.7 


2.0 


2.2 


2.6 


2.8 



O- 




-6 -4-2 2 

ENERGY (eV) 



-4-2 2 

ENERGY (eV) 



FIG. 3. (Color online) Comparison of the partial spectral 
functions obtained from the LDA+t/ (gray lines) and Ander- 
son model calculations for the Mn/CuN (left) and Fe/CuN 
(right) systems {U=3 eV and Jh= 0.9 eV). 
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FIG. 4. (Color online) Comparison of the partial spectral 
functions obtained from the LDA+[/ (gray lines) and Ander- 
son model calculations for the Co/CuN system ([/=3 eV and 
Jh= 0.9 eV). 



response theory, 
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where En and \n) are the eigenvalues and eigenstates of 
the Anderson Hamiltonian, and A^'n = En' — En is an 
excitation energy. The corresponding spectral functions 
are antisymmetric^ i.e. they are positive for cj > and 
negative for cj < 0^ Here we phenomenologically intro- 
duce the finite width for the collective spin excitations, 
r(A^^/) that relates to the lifetimes of these excitations, 
r = ^ . This allows us to decrease the effects of the dis- 
cretized bath on the excitation spectrum. In general, the 
width can be calculated by using the self-energy part of 
the two-particle Green's function. To simplify the con- 
sideration, we use the principal property of T, that is 



r(Aw) « A^ 



(4) 



The corresponding estimation for the ratio ^^ "^"^'^^ can 
be found in the previous theoretical works. For instance, 
the authors of Refo have fitted the experimental STM 
curve for Fe/CuN with the width F = 0.01 - 0.03 meV. 
Taking into account the spin excitation energy of 1 meV, 



we can estimate 



r(A,./). 



0.01. 



Since we simulate the magnetic states of the impurity 
with magnetic moments ranging from 2.2 eV (Co/CuN) 



to 4.4 eV (Mn/CuN), the intensity of the z2;-component 
of the response function is strongly suppressed. Below 
we will analyze the xf ^^^ Xi contributions that are 
related to each other by Im{Xi{~^)) ~ ~-^'^{xf{'^))' 

The resulting spectral functions are presented in Fig. 5. 
One can see that there are well-separated low- and 
high-energy excitations with the gap ranging from 3 eV 
(Mn/CuN) to 1 eV (Co/CuN). The low-energy peak is 
a response of the system on the external magnetic field, 
Beff, as the resonance energy is fully controlled by the 
static 2:-component of the magnetic field, Bq. At the 
same time, the high-energy structure corresponds to the 
intra-atomic spin-flip excitations. Analyzing these exci- 
tations at the time scale, it can be seen the intra-atomic 
excitations (1-2 eV) are typical for ultrafast processes at 
the femtosecond scale, whereas the response to the exter- 
nal magnetic field can be associated with the nanoscale 
range. 

Importantly, all the magnetic orbitals of Fe and Mn 
atoms demonstrate the same intensity of the low-energy 
peak. It is not the case for Co impurity where the inten- 
sity of the response function depends on the symmetry of 
the 3d orbital with the largest response resulting from the 
x'^ — y'^ orbital. The smallest intensity of the low-energy 
peak is demonstrated by the xy orbital. The difference 
between them is due to the hybridization with the nitro- 
gen states. While the orbital of the xy symmetry points 
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FIG. 5. (Color online) Low- (left panel) and high-energy 
(right panel) parts of the spin-flip response functions (in 
{giiiB)'^eV~^) calculated with [7=3 eV and Jh= 0.9 eV. The 
blue, red, purple and brown lines correspond to the xy, yz 
(xz), 3z^ — r^ and x^ — y^ orbitals. The gray lines are exam- 
ples of the fitting by using Bloch equation, Eq.®. 



TABLE III. Orbital relaxation times, T2 estimated by mini- 
mizing Eq.(6). 



xy yz 


3z'^ -r 


^ xz x^ — y^ 


Low-energy 


(slow) 


excitations 


Mn 7 ps 7 ps 
Fe 6 ps 7 ps 
Co 40 fs 50 fs 


7 ps 
7 ps 
60 fs 


7 ps 7 ps 
7 ps 7 ps 
50 fs 80 fs 


High-energ}^ 


(fast) 


excitations 


Mn 70 fs 70 fs 
Fe 17 fs 22 fs 
Co 9 fs 11 fs 


70 fs 
15 fs 
9fs 


70 fs 70 fs 
22 fs 11 fs 
11 fs 12 fs 



in the direction of the largest overlap with nitrogen, the 
x^ — y^ orbital lobes look away from nitrogen. 

Orbital relaxation times. For practical purposes it is 
important to show a connection between the lifetimes of 
quantum excitations and relaxation times of the mag- 
netic moments of individual impurity orbitals. For this 
purpose, we use the Bloch's^^ approach and consider 
the nonequilibrium behavior of the orbital magnetic mo- 
ments of the impurity in the magnetic field, Beff. This 
field has a stationary z-component of frequency ujq = 
qIIbBq and weak time-dependent component along x 
axis. Thus B^ff = (^i cos uot^ 0, ^o)- 

The impurity in the strong magnetic field applied along 
z direction yields a finite magnetization. The z compo- 
nent of the magnetization of the \th orbital, ml fluctu- 
ates near stationary solution m^. Taking into account a 
small oscillating magnetic field in the x-direction we ob- 
tain that the x and y components of the magnetization 
oscillate near the zero value. If mf is perturbed from the 
equilibrium state {mf^m^ ^ 0) the return to this state 
will take characteristic times T^ (longitudinal) and T2 
(transverse). The corresponding equations of motion for 
the X and y components of the impurity orbital magne- 
tization are given by 



dm] 



x,y 



dt 



[m X Beff]cc,y 



(5) 



The solution of this system of differential equations cor- 
responds to the following susceptibility— 
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where m!^ is the equilibrium value of the magnetic mo- 
ment of the \th orbital. Approximating nonequilibrium 
Xf^(co') by the spin- flip susceptibilities obtained within 
the linear response theory, xf^(^) = i(Xi^ +Xi ) we can 
estimate the orbital relaxation times. 

The spin-flip response function contains low- and high- 
energy structures that correspond to different types of 
spin excitation in the system. We substitute these parts 
into Eq.(6) separately to simulate the magnetization dy- 
namics on different times scales. In this case, the pa- 
rameter m!l controls the capacity of the corresponding 
susceptibility part and it can only be associated with 




FIG. 6. (Color online) Schematic representation of the Mn 
dimer on the CuN surface (top view) simulated in this work. 



the equilibrium magnetic moment of the impurity orbital. 
Let us start from the analysis of the low-energy excita- 
tions attributed to a response to the external magnetic 
field. The obtained relaxation times are presented in Ta- 
ble III. For Mn/CuN and Fe/CuN systems, we obtain the 
same relaxation times for the magnetic orbit als. In case 
of cobalt system, there is a small variation of T2 for dif- 
ferent impurity orbital states. Another important result 
is that the cobalt relaxation times are of two orders of 
magnitude faster than those of iron and manganese sys- 
tems. It is due to the fact that the Co/CuN system is 
more itinerant than the others. 

The spin-flip excitations of intra-atomic nature corre- 
spond to the high-energy part of the spectral function for 
uj > 1 eV. The minimization procedure for this part of the 
spectrum leads to a picture of ultrafast relaxation times 
for the orbital magnetization (Table IV). The main differ- 
ence from the low-energy excitation is that the Fe/CuN 
and Co/CuN systems demonstrate a complex cascade of 
the spin-flip excitations within the energy window of 1-2 
eV. The manganese impurity shows two isolated excita- 
tions. 

Correlated exchange in Mn- dimer. The next step of 
our study is to examine the influence of the dynamical 
correlations on the isotropic exchange interactions be- 
tween adatoms on the surface. For these purposes we 
chose the Mn-dimer on the CuN surface^^ that demon- 
strates the largest exchange interaction of 6.4 meV. The 
geometry of the Mn-dimer on the CuN surface is de- 
scribed in detail in our previous investigation.^^ To define 
the parameters of the two-impurity Anderson Hamilto- 
nian, we use a combination of the projection procedure 
and minimization of the on-site LDA Green's function. 
The detailed discussion is presented in the Supplemen- 
tary materials. The constructed model is solved by us- 
ing an exact diagonalization technique. The comparison 
of the calculated spectral functions with those obtained 
from the LDA-\-U calculations is presented in Fig. 7. One 
can see that, in contrast to LDA-\-U results, the Anderson 
model demonstrates additional excitations at the Fermi 
level for the 3z^ — r^ and xz orbitals. 

Before our discussing the two-impurity Anderson 
model and the results of its solution, let us analyze quali- 
tatively possible effects of the dynamical Coulomb corre- 
lations on the exchange interaction. For that we consider 
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FIG. 7. (Color online) Left panel. Comparison of the par- 
tial spectral functions obtained from LDA (gray lines) and 
non-interacting model Hamiltonian for the Mn dimer on CuN 
system. Right panel. Comparison of the partial spectral func- 
tions obtained from LDA+?7 (gray lines) and Anderson model 
calculations for the Mn dimer on the CuN surface with t/=3 
eV and Jh^ 0.9 eV. 



a general expression for the total exchange interaction of 



a given site with all magnetic environment ^i2ii2S 

1 f^^ 

— /m / AoGioAoG^oduj, (7) 



where Gao is the spin-dependent on-site local Green's 
function, the on-site potential Aq = G7q — GIq and Ef 
is the Fermi energy. The Green's function taking into 
account the dynamical Coulomb correlations that can 
be obtained, for instance, in the LSDA+S approach^ 
is written in the following form: 



Gao(^) 



k 



(a; + M-e'^(k)-S^H)- 



(8) 



where a denotes spin, /i is the chemical potential, Sq is 
the self-energy and e^(k) is the LSDA spectrum. If the 
self-energy part in Eq.(j8]) is zero then we obtain the local 
LSDA Green's function. 

Let us compare the magnetic couplings calculated 
in LSDA (static mean-field) and LSDA+S (dynamical 
mean- field) approaches. For that, we expand the real 
part of the self-energy in the vicinity of the Fermi energy 
leaving only linear term, we obtain 



where the quasiparticle residue is given by 

dRe^juj) , 1 



(9) 



(10) 



In turn, the correlated Green's function, Eq.(j8]) can be 
reduced to the effective LSDA form 



G.oH«E(c^^-'-e'^(k))-\ 



(11) 



of the magnetic coupling within the Anderson model ap- 
proach, we used the exchange interaction in a strongly- 
correlated formulation^ 



Ji2 = ^Irnj " ^IGl^^lGi.doj, 



(13) 



where S^ = Sj — T>j is the self-energy part of the on- 
site Green's function and G12 is the inter-site correlated 
Green's function. 

We now turn to the analysis of the calculated exchange 
interaction between manganese impurities on CuN. In the 
previous work^ we demonstrated the effect of the static 
Coulomb correlations by comparing magnetic couplings 
calculated by employing LSDA and LDA+/7 approaches. 
The account of the U correction leads to a strong renor- 
malization of the exchange interaction. Here one can see 
that the account of the dynamical correlations leads to 
an additional decrease of the exchange interaction be- 
tween magnetic moments of the manganese atoms. We 
note that the approximately the same value of the to- 
tal exchange interaction is obtained in the LDA+[/ and 
Anderson model calculations with completely different 
values of the on-site Coulomb interaction, U=3 eV (An- 
derson) and [7=5.88 eV (LDA+U).^^ 
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with rescaled frequencies. This correlated Green's func- 
tion being substituted into Eq.Q leads us to the follow- 
ing ratio 



jLSDA^J: _ ry jLSDA 
Jq — ^-'0 



(12) 



The performed qualitative consideration demonstrates 
that many-body correlation effects tend to renormalize a 
static mean-field (LSDA) value of the isotropic magnetic 
coupHng. It is easy to show that the anisotropic inter- 
actions, such as Dzyaloshinskii-Moriya interaction and 
single-ion anisotropy are satisfied to the same relation 
(p!2|) between static and dynamical values. Quantities 
related with the exchange interaction are also rescaled. 
For instance, one can consider the Curie- Weiss temper- 
ature that in the high-temperature limit can be written 
as Tew 



JoSjS+l) 
SkE 



Since the impurity systems we study demonstrate a 
complex many-orbital behavior, it requires accurate nu- 
merical calculations. For reliable quantitative estimates 



FIG. 8. (Color online) Exchange interactions calculated in 
LDA+[/ ([7=5.88 eV) and Anderson model ([7 = 3 eV) as a 
function of the Fermi energy. 

The decrease of the orbital contributions to the ex- 
change interaction in LDA+/7 and Anderson model cal- 
culations with respect to the LSDA results is not uni- 
form. As can be seen from Table IV the interaction be- 
tween the xy orbitals of manganese atoms demonstrates 
the strongest renormalization. This allows us to conclude 
that the account of the frequency dependence of the self- 
energy plays a crucial role for the correct description of 
the magnetic couplings in surface nanostructures. The 
best agreement with the experimental value of 6.4 meV 
is achieved by using [7=3.1 eV in the Anderson model 
calculations. 

Since applying the voltage between the tip and sam- 
ple in STM experiments one shifts the Fermi level, Ep 
of the sample, it is of interest to examine the stability of 
the antiferromagnetic coupling with respect to Ep. Fig. 8 



TABLE IV. Total exchange interaction, J12 and the largest 
orbital contributions (in meV) to the magnetic coupling ob- 
tained from the LSDA, LDA+U ([/=5.88 eV)^^ and Anderson 
model ([/=3 eV) calculations. 



Clxy~Clxy ^3^2 



Jl2 



LSDA 
LDA+U 
Anderson 



5.9 
L8 
L5 



L7 
L2 
L4 



2.7 
1.7 
1.1 



15.3 
6.5 
7.1 



gives the dependence of the exchange interaction on the 
upper integration limit in Eq. (p!3|) . One can see that the 
deviation from the half-filled state in Mn-dimer leads to 
ferromagnetic solutions both in the LDA+/7 and Ander- 
son model calculations. 

To conclude, using the combination of the LDA and 
Anderson model calculations we have studied the elec- 
tronic and magnetic excitations in Mn, Fe, and Co im- 
purities deposited on the CuN surface. In the frame- 
work of the linear response theory, we observe fast intra- 
atomic spin-flip excitations at high energies and slow ex- 



citation that are due to the external magnetic field. The 
estimated relaxation times of the magnetic moments of 
3d- Co orbit als are found to be two times smaller than 
that for iron and manganese adatoms on CuN. The two- 
impurity model calculations reveal a strong renormaliza- 
tion of of the exchange interaction between magnetic mo- 
ments of Mn atoms due to the dynamical correlation ef- 
fects. This result should be carefully taken into account 
for a realistic modeling of surface nanostructures. 
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